Sequential estimation of neural models by Bayesian filtering

Un dels reptes més difícils de la neurociència és el d'entendre la connectivitat del cervell. Aquest problema es pot tractar des de diverses perspectives, aquí ens centrem en els fenòmens locals que ocorren en una sola neurona. L'objectiu final és, doncs, entendre la dinàmica de les neurones i com la interconnexió amb altres neurones afecta al seu estat. Les observacions de traces del potencial de membrana constitueixen la principal font d'informació per a derivar models matemàtics d'una neurona, amb cert sentit biofísic. En particular, la dinàmica de les variables auxiliars i els paràmetres del model són estimats a partir d'aquestes traces de voltatge. El procés és en general costós i típicament implica una gran varietat de blocatges químics de canals iònics, així com una certa incertesa en els valors dels paràmetres a causa del soroll de mesura. D'altra banda, les traces de potencial de membrana també són útils per obtenir informació valuosa sobre l'entrada sinàptica, un problema invers sense solució satisfactòria a hores d'ara. En aquesta Tesi, estem interessats en mètodes d'estimació seqüencial, que permetin evitar la necessitat de repeticions que podrien ser contaminades per la variabilitat neuronal. En particular, ens concentrem en mètodes per extreure l'activitat intrínseca dels canals iònics, és a dir, les probabilitats d'obertura i tancament de canals iònics, i la contribució de les conductàncies sinàptiques. Hem dissenyat un mètode basat en la teoria Bayesiana de filtrat per inferir seqüencialment aquestes quantitats a partir d'una única traça de voltatge, potencialment sorollosa. El mètode d'estimació proposat està basat en la suposició d'un model de neurona conegut. Això és cert fins a cert punt, però la majoria dels paràmetres en el model han de ser estimats per endavant (això és valid per a qualsevol model). Per tant, el mètode s'ha millorat pel cas de models amb paràmetres desconeguts, incloent-hi un procediment per estimar conjuntament els paràmetres i les variables dinàmiques. Hem validat els mètodes d'inferència proposats mitjançant simulacions realistes. Les prestacions en termes d'error d'estimació s'han comparat amb el límit teòric, que s'ha derivat també en el marc d'aquesta Tesi

Maximum likelihood estimation of position in GNSS

In this letter, we obtain the Maximum Likelihood Estimator of position in the framework of Global Navigation Satellite Systems. This theoretical result is the basis of a completely different approach to the positioning problem, in contrast to the conventional two-steps position estimation, consisting of estimating the synchronization parameters of the in-view satellites and then performing a position estimation with that information. To the authors’ knowledge, this is a novel approach which copes with signal fading and it mitigates multipath and jamming interferences. Besides, the concept of Position–based Synchronization is introduced, which states that synchronization parameters can be recovered from a user position estimation. We provide computer simulation results showing the robustness of the proposed approach in fading multipath channels. The Root Mean Square Error performance of the proposed algorithm is compared to those achieved with state-of-the-art synchronization techniques. A Sequential Monte–Carlo based method is used to deal with the multivariate optimization problem resulting from the ML solution in an iterative way.Peer Reviewe

A Game Theoretical algorithm for joint power and topology control in distributed WSN

In this paper, the issue of network topology control in wireless networks using a fully distributed algorithm is considered. Whereas the proposed distributed algorithm is designed applying game theory concepts to design a non-cooperative game, network connectivity is guaranteed based on asymptotic results of network connectivity. Simulations show that for a relatively low node density, the probability that the proposed algorithm leads to a connected network is close to one.Peer ReviewedPostprint (published version

Benchmarking CPUs and GPUs on embedded platforms for software receiver usage

Smartphones containing multi-core central processing units (CPUs) and powerful many-core graphics processing units (GPUs) bring supercomputing technology into your pocket (or into our embedded devices). This can be exploited to produce power-efficient, customized receivers with flexible correlation schemes and more advanced positioning techniques. For example, promising techniques such as the Direct Position Estimation paradigm or usage of tracking solutions based on particle filtering, seem to be very appealing in challenging environments but are likewise computationally quite demanding. This article sheds some light onto recent embedded processor developments, benchmarks Fast Fourier Transform (FFT) and correlation algorithms on representative embedded platforms and relates the results to the use in GNSS software radios. The use of embedded CPUs for signal tracking seems to be straight forward, but more research is required to fully achieve the nominal peak performance of an embedded GPU for FFT computation. Also the electrical power consumption is measured in certain load levels.Peer ReviewedPostprint (published version

Feasibility of snapshot GNSS carrier phase amendment based on LAMBDA ratio tests

This research was funded by Albora Technologies and Universitat Politècnica de Catalunya with industrial PhD grant number DI 082 from the Generalitat de Catalunya; This research was partially funded by the Spanish Ministry of Science and Innovation project RTI2018-094295-B-I00. P.C. has been partially supported by the NSF under Awards CNS-1815349 and ECCS-1845833.Peer ReviewedPostprint (author's final draft

A Method to determine secondary codes and carrier phases of short snapshot signals

Recently, the Snapshot Real-Time Kinematic (SRTK) technique was demonstrated, which aims at achieving high accuracy navigation solutions with a very short signal collection. The main challenge in implementing SRTK is the generation of valid carrier-phase measurements, which relies on a data bit ambiguity (DBA) resolution process. For pilot signals, this step is equivalent to the correct selection of secondary code indexes (SCIs) from the ambiguous sets obtained from a multi-hypotheses (MH) acquisition process. Currently, SCI ambiguities are solved independently for each satellite. However, this method is ineffective when the snapshot signal is relatively short. In order to tackle this problem, this article proposes a new method that makes use of assistance data and processes information from all satellites to jointly solve the DBA issue. This new method is shown to be more effective in determining the correct SCI and enabling valid snapshot carrier-phase measurements, largely expanding the scope of high-accuracy snapshot positioning.This research was supported by the Albora Technologies and Universitat Politècnica de Catalunya with industrial PhD grant number DI 082 from the Generalitat de Catalunya and the project RTI2018-094295-B-I00 funded by the MCIN/AEI 10.13039/501100011033 which is co-funded by the FEDER programme. P.C. has been partially supported by the NSF under Awards CNS-1815349 and ECCS-1845833.Peer ReviewedPostprint (published version

Cloud-based single-frequency Snapshot RTK positioning

With great potential for being applied to Internet of Things (IoT) applications, the concept of cloud-based Snapshot Real Time Kinematics (SRTK) was proposed and its feasibility under zero-baseline configuration was confirmed recently by the authors. This article first introduces the general workflow of the SRTK engine, as well as a discussion on the challenges of achieving an SRTK fix using actual snapshot data. This work also describes a novel solution to ensure a nanosecond level absolute timing accuracy in order to compute highly precise satellite coordinates, which is required for SRTK. Parameters such as signal bandwidth, integration time and baseline distances have an impact on the SRTK performance. To characterize this impact, different combinations of these settings are analyzed through experimental tests. The results show that the use of higher signal bandwidths and longer integration times result in higher SRTK fix rates, while the more significant impact on the performance comes from the baseline distance. The results also show that the SRTK fix rate can reach more than 93% by using snapshots with a data size as small as 255 kB. The positioning accuracy is at centimeter level when phase ambiguities are resolved at a baseline distance less or equal to 15 km.This research was funded by Albora Technologies and Universitat Politècnica de Catalunya with industrial PhD grant number DI 082 from the Generalitat de Catalunya; This research was partially funded by the Spanish Ministry of Science and Innovation project RTI2018-094295-B-I00. P.C. has been partially supported by the NSF under Awards CNS-1815349 and ECCS-1845833Peer ReviewedPostprint (published version

Bayesian signal processing techniques for GNSS receivers: from multipath mitigation to positioning

Aquesta tesi gira al voltant del disseny de receptors per a sistemes globals de navegació per satèl·lit (Global Navigation Satellite Systems, GNSS). El terme GNSS fa referència a tots aquells sistemes de navegació basats en una constel·lació de satèl·lits que emeten senyals de navegació útils per a posicionament. El més popular és l'americà GPS, emprat globalment. Els esforços d'Europa per a tenir un sistema similar veuran el seu fruit en un futur proper, el sistema s'anomena Galileo. Altres sistemes globals i regionals existeixen dissenyats per al mateix objectiu: calcular la posició dels receptors. Inicialment la tesi presenta l'estat de l'art en GNSS, a nivell de l'estructura dels actuals senyals de navegació i pel que fa a l'arquitectura dels receptors.El disseny d'un receptor per a GNSS consta d'un seguit de blocs funcionals. Començant per l'antena receptora fins al càlcul final de la posició del receptor, el disseny proporciona una gran motivació per a la recerca en diversos àmbits. Tot i que la cadena de Radiofreqüència del receptor també és comentada a la tesis, l'objectiu principal de la recerca realitzada recau en els algorismes de processament de senyal emprats un cop realitzada la digitalització del senyal rebut. En un receptor per a GNSS, aquests algorismes es poden dividir en dues classes: els de sincronisme i els de posicionament. Aquesta classificació correspon als dos grans processos que típicament realitza el receptor. Primer, s'estima la distancia relativa entre el receptor i el conjunt de satèl·lits visibles. Aquestes distancies es calculen estimant el retard patit pel senyal des de que és emès pel corresponent satèl·lit fins que és rebut pel receptor. De l'estimació i seguiment del retard se n'encarrega l'algorisme de sincronisme. Un cop calculades la distancies relatives als satèl·lits, multiplicant per la velocitat de la llum el retards estimats, l'algorisme de posicionament pot operar. El posicionament es realitza típicament pel procés de trilateralització: intersecció del conjunt d'esferes centrades als satèl·lits visibles i de radi les distancies estimades relatives al receptor GNSS. Així doncs, sincronització i posicionament es realitzen de forma seqüencial i ininterrompudament. La tesi fa contribucions a ambdues parts, com explicita el subtítol del document.Per una banda, la tesi investiga l'ús del filtrat Bayesià en el seguiment dels paràmetres de sincronisme (retards, desviaments Doppler i phases de portadora) del senyal rebut. Una de les fonts de degradació de la precisió en receptors GNSS és la presència de repliques del senyal directe, degudes a rebots en obstacles propers. És per això que els algorismes proposats en aquesta part de la tesi tenen com a objectiu la mitigació de l'efecte multicamí. La dissertació realitza una introducció dels fonaments teòrics del filtrat Bayesià, incloent un recull dels algorismes més populars. En particular, el Filtrat de Partícules (Particle Filter, PF) s'estudia com una de les alternatives més interessants actualment per a enfrontar-se a sistemes no-lineals i/o no-Gaussians. Els PF són mètodes basats en el mètode de Monte Carlo que realitzen una caracterització discreta de la funció de probabilitat a posteriori del sistema. Al contrari d'altres mètodes basats en simulacions, els PF tenen resultats de convergència que els fan especialment atractius en casos on la solució òptima no es pot trobar. En aquest sentit es proposa un PF que incorpora un seguit de característiques que el fan assolir millors prestacions i robustesa que altres algorismes, amb un nombre de partícules reduït. Per una banda, es fa un seguiment dels estats lineals del sistema mitjançant un Filtre de Kalman (KF), procediment conegut com a Rao-Blackwellization. Aquest fet provoca que la variància de les partícules decreixi i que un menor nombre d'elles siguin necessàries per a assolir una certa precisió en l'estimació de la distribució a posteriori. D'altra banda, un dels punts crítics en el disseny de PF és el disseny d'una funció d'importància (emprada per a generar les partícules) similar a l'òptima, que resulta ésser el posterior. Aquesta funció òptima no està disponible en general. En aquesta tesi, es proposa una aproximació de la funció d'importància òptima basada en el mètode de Laplace. Paral·lelament es proposen algorismes com l'Extended Kalman Filter (EKF) i l'Unscented Kalman Filter (UKF), comparant-los amb el PF proposat mitjançant simulacions numèriques.Per altra banda, la presentació d'un nou enfocament al problema del posicionament és una de les aportacions originals de la tesi. Si habitualment els receptors operen en dos passos (sincronització i posicionament), la proposta de la tesi rau en l'Estimació Directa de la Posició (Direct Position Estimation, DPE) a partir del senyal digital. Tenint en compte la novetat del mètode, es proporcionen motivacions qualitatives i quantitatives per a l'ús de DPE enfront al mètode convencional de posicionament. Se n'ha estudiat l'estimador de màxima versemblança (Maximum Likelihood, ML) i un algorisme per a la seva implementació pràctica basat en l'algorisme Accelerated Random Search (ARS). Els resultats de les simulacions numèriques mostren la robustesa de DPE a escenaris on el mètode convencional es veu degradat, com per exemple el cas d'escenaris rics en multicamí. Una de les reflexions fruit dels resultats és que l'ús conjunt dels senyals provinents dels satèl·lits visibles proporciona millores en l'estimació de la posició, doncs cada senyal està afectada per un canal de propagació independent. La tesi també presenta l'extensió de DPE dins el marc Bayesià: Bayesian DPE (BDPE). BDPE manté la filosofia de DPE, tot incloent-hi possibles fonts d'informació a priori referents al moviment del receptor. Es comenten algunes de les opcions com l'ús de sistemes de navegació inercials o la inclusió d'informació atmosfèrica. Tot i així, cal tenir en compte que la llista només està limitada per la imaginació i l'aplicació concreta on el marc BDPE s'implementi.Finalment, la tesi els límits teòrics en la precisió dels receptors GNSS. Alguns d'aquests límits teòrics eren ja coneguts, d'altres veuen ara la llum. El límit de Cramér-Rao (Cramér-Rao Bound, CRB) ens prediu la mínima variància que es pot obtenir en estimar un paràmetre mitjançant un estimador no esbiaixat. La tesi recorda el CRB dels paràmetres de sincronisme, resultat ja conegut. Una de les aportacions és la derivació del CRB de l'estimador de la posició pel cas convencional i seguint la metodologia DPE. Aquests resultats proporcionen una comparativa asimptòtica dels dos procediments pel posicionament de receptors GNSS. D'aquesta manera, el CRB de sincronisme pel cas Bayesià (Posterior Cramér-Rao Bound, PCRB) es presenta, com a límit teòric dels filtres Bayesians proposats en la tesi.This dissertation deals with the design of satellite-based navigation receivers. The term Global Navigation Satellite Systems (GNSS) refers to those navigation systems based on a constellation of satellites, which emit ranging signals useful for positioning. Although the american GPS is probably the most popular, the european contribution (Galileo) will be operative soon. Other global and regional systems exist, all with the same objective: aid user's positioning. Initially, the thesis provides the state-of-the-art in GNSS: navigation signals structure and receiver architecture. The design of a GNSS receiver consists of a number of functional blocks. From the antenna to the final position calculation, the design poses challenges in many research areas. Although the Radio Frequency chain of the receiver is commented in the thesis, the main objective of the dissertation is on the signal processing algorithms applied after signal digitation. These algorithms can be divided into two: synchronization and positioning. This classification corresponds to the two main processes typically performed by a GNSS receiver. First, the relative distance between the receiver and the set of visible satellites is estimated. These distances are calculated after estimating the delay suffered by the signal traveling from its emission at the corresponding satellite to its reception at the receiver's antenna. Estimation and tracking of these parameters is performed by the synchronization algorithm. After the relative distances to the satellites are estimated, the positioning algorithm starts its operation. Positioning is typically performed by a process referred to as trilateration: intersection of a set of spheres centered at the visible satellites and with radii the corresponding relative distances. Therefore, synchronization and positioning are processes performed sequentially and in parallel. The thesis contributes to both topics, as expressed by the subtitle of the dissertation.On the one hand, the thesis delves into the use of Bayesian filtering for the tracking of synchronization parameters (time-delays, Doppler-shifts and carrier-phases) of the received signal. One of the main sources of error in high precision GNSS receivers is the presence of multipath replicas apart from the line-of-sight signal (LOSS). Wherefore the algorithms proposed in this part of the thesis aim at mitigating the multipath effect on synchronization estimates. The dissertation provides an introduction to the basics of Bayesian filtering, including a compendium of the most popular algorithms. Particularly, Particle Filters (PF) are studied as one of the promising alternatives to deal with nonlinear/nonGaussian systems. PF are a set of simulation-based algorithms, based on Monte-Carlo methods. PF provide a discrete characterization of the posterior distribution of the system. Conversely to other simulation-based methods, PF are supported by convergence results which make them attractive in cases where the optimal solution cannot be analytically found. In that vein, a PF that incorporates a set of features to enhance its performance and robustness with a reduced number of particles is proposed. First, the linear part of the system is optimally handled by a Kalman Filter (KF), procedure referred to as Rao-Blackwellization. The latter causes a reduction on the variance of the particles and, thus, a reduction on the number of required particles to attain a given accuracy when characterizing the posterior distribution. A second feature is the design of an importance density function (from which particles are generated) close to the optimal, not available in general. The selection of this function is typically a key issue in PF designs. The dissertation proposes an approximation of the optimal importance function using Laplace's method. In parallel, Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF) algorithms are considered, comparing these algorithms with the proposed PF by computer simulations.On the other hand, a novel point of view in the positioning problem constitutes one of the original contributions of the thesis. Whereas conventional receivers operate in a two-steps procedure (synchronization and positioning), the proposal of the thesis is a Direct Position Estimation (DPE) from the digitized signal. Considering the novelty of the approach, the dissertation provides both qualitative and quantitative motivations for the use of DPE instead of the conventional two-steps approach. DPE is studied following the Maximum Likelihood (ML) principle and an algorithm based on the Accelerated Random Search (ARS) is considered for a practical implementation of the derived estimator. Computer simulation results carried show the robustness of DPE in scenarios where the conventional approach fails, for instance in multipath-rich scenarios. One of the conclusions of the thesis is that joint processing of satellite's signals provides enhance positioning performances, due to the independent propagation channels between satellite links. The dissertation also presents the extension of DPE to the Bayesian framework: Bayesian DPE (BDPE). BDPE maintains DPE's philosophy, including the possibility of accounting for sources of side/prior information. Some examples are given, such as the use of Inertial Measurement Systems and atmospheric models. Nevertheless, we have to keep in mind that the list is only limited by imagination and the particular applications were BDPE is implemented. Finally, the dissertation studied the theoretical lower bounds of accuracy of GNSS receivers. Some of those limits were already known, others see the light as a result of the research reported in the dissertation. The Cramér-Rao Bound (CRB) is the theoretical lower bound of accuracy of any unbiased estimator of a parameter. The dissertation recalls the CRB of synchronization parameters, result already known. A novel contribution ofthe thesis is the derivation of the CRB of the position estimator for either conventional and DPE approaches. These results provide an asymptotical comparison of both GNSS positioning approaches. Similarly, the CRB of synchronization parameters for the Bayesian case (Posterior Cramér-Rao Bound, PCRB) is given, used as a fundamental limit of the Bayesian filters proposed in the thesis

Sequential estimation of neural models by Bayesian filtering

Un dels reptes més difícils de la neurociència és el d'entendre la connectivitat del cervell. Aquest problema es pot tractar des de diverses perspectives, aquí ens centrem en els fenòmens locals que ocorren en una sola neurona. L'objectiu final és, doncs, entendre la dinàmica de les neurones i com la interconnexió amb altres neurones afecta al seu estat. Les observacions de traces del potencial de membrana constitueixen la principal font d'informació per a derivar models matemàtics d'una neurona, amb cert sentit biofísic. En particular, la dinàmica de les variables auxiliars i els paràmetres del model són estimats a partir d'aquestes traces de voltatge. El procés és en general costós i típicament implica una gran varietat de blocatges químics de canals iònics, així com una certa incertesa en els valors dels paràmetres a causa del soroll de mesura. D'altra banda, les traces de potencial de membrana també són útils per obtenir informació valuosa sobre l'entrada sinàptica, un problema invers sense solució satisfactòria a hores d'ara. En aquesta Tesi, estem interessats en mètodes d'estimació seqüencial, que permetin evitar la necessitat de repeticions que podrien ser contaminades per la variabilitat neuronal. En particular, ens concentrem en mètodes per extreure l'activitat intrínseca dels canals iònics, és a dir, les probabilitats d'obertura i tancament de canals iònics, i la contribució de les conductàncies sinàptiques. Hem dissenyat un mètode basat en la teoria Bayesiana de filtrat per inferir seqüencialment aquestes quantitats a partir d'una única traça de voltatge, potencialment sorollosa. El mètode d'estimació proposat està basat en la suposició d'un model de neurona conegut. Això és cert fins a cert punt, però la majoria dels paràmetres en el model han de ser estimats per endavant (això és valid per a qualsevol model). Per tant, el mètode s'ha millorat pel cas de models amb paràmetres desconeguts, incloent-hi un procediment per estimar conjuntament els paràmetres i les variables dinàmiques. Hem validat els mètodes d'inferència proposats mitjançant simulacions realistes. Les prestacions en termes d'error d'estimació s'han comparat amb el límit teòric, que s'ha derivat també en el marc d'aquesta Tesi

Electrònica de les radiocomunicacions navals

2017/201